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Введение. Рассматривается новый вид устройств для сбора 
и аккумулирования энергии — воздушный надувной 
коллектор. Как 
устанавливаются стационарно, что не подразумевает 


правило, надувные — коллекторы 
ориентацию коллектора вслед за движением Солнца. Ввиду 


НИЗКОЙ себестоимости — рассматриваемых = изделий, 
необходимо предложить и исследовать максимально 
эффективную конструкцию. 
Материалы и методы. Рассматривается коллектор, 
состоящий из последовательно соединенных замкнутых 
цилиндрических сегментов. Полости цилиндров заполнятся 
воздухом, что обеспечивает постоянство конструкции. 
Математическое моделирование определения 
температурного поля воздушного надувного коллектора 
выполнено с помощью метода конечных элементов. 
Результаты исследования. Распределение температурного 
поля, в зависимости от направления потока солнечной 
радиации, были подтверждены = экспериментально. 
Математические модели признаны адекватными. Прирост 
температуры теплопоглощающего слоя по отношению к 
температуре окружающей среды составил от 7° до 26,2° в 
зависимости от части сегмента коллектора. 

Обсуждение и заключения. В ходе имитационных и 
экспериментальных исследований было установлено, что на 
эффективность коллектора, кроме инсоляции Солнца, 
оказывают воздействие комплекс факторов окружающей 
среды: влажность воздуха, сила ветра и др. При этом 
необходимо исполнение 


устройства, а также эксплуатационные и теплофизические 


учитывать конструктивное 


характеристики Детальное 


исследование степени воздействия внешних и внутренних 


применяемых материалов. 


факторов на температурное поле коллектора требует 
последующей разработки программного комплекса. 
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Введение. Популяризация применения альтернативных источников энергии требует детального 
изучения новых гелиотехнических решений. Тем более, если речь идет о малоизученном сегменте надувных 
воздушных коллекторов (далее коллектор). Сложность и многофакторность процессов теплообмена в 
рассматриваемых системах определяют трудобёмкость сбора и анализа экспериментальных данных [1-8]. На 
этапе разработки новых устройств, а также оптимизации действующей модели коллектора, важна аналитика по 
каждому фактору, влияющему на эффективность устройства в целом. Разработанная математическая модель 
позволит изучить температурное поле коллектора в зависимости от направления потока солнечной радиации. 

Постановка задачи. В качестве объекта исследования был выбран коллектор, представляющий собой 
систему последовательно соединенных замкнутых цилиндров, наполненных воздухом. Каждый цилиндр 
системы состоит из двух слоев: верхнего — светопропускающего и нижнего — теплопоглощающего 
(абсорбера). Касательно к выступающим точкам теплопоглощающей части цилиндра прикреплена 
теплоизоляционная подложка. В образовавшемся зазоре локализуется теплоноситель — воздух. Предлагаемая 
конструкция коллектора (рис. 1) позволяет обеспечить стабильно высокую температуру теплоносителя вне 
зависимости от направления потока солнечной радиации на поверхность коллектора. 

Результаты натурных и лабораторных экспериментов [9] определили основные требования, 
предъявляемые к материалам коллектора. Для изготовления светопропускающего слоя коллектора был выбран 
полиэтилентерефталат. Данный материал обладает высокими показателями теплостойкости, стойкости к УФ- 
старению, а также коэффициентом светопропускаемости равным 0,85. Теплопоглощающий слой выполнен из 
шерстяного полотна, оно вкладывается по контуру нижнего полуцилиндра. 

Рассмотрим процесс теплообмена в одном звене солнечного коллектора (рис. 2) в зависимости от 
направления потока радиации на поверхность коллектора. Ввиду наличия прозрачного верхнего слоя 
предложено допущение, что верхний слой не оказывает значительного воздействия на распределение 
температурного поля цилиндрического сегмента. 





Рис. 1. Изображение воздушного надувного коллектора в разрезе: 
| — цилиндрический сегмент коллектора; 2 — светопропускающий прозрачный слой сегмента; 


3 — теплопоглощающий слой сегмента; 
4 — теплоизоляционный слой, 5 — область локализации теплоносителя 


Ето. 1. РтоШе ореаг-тПаюЫе соПес1ог: 
1 - суппаяса[ соПесог 5езтеии; 2 - тапуисеп! тапзрагеи! зестеп! [ауег; 
3 - пеа-аб5огЬте 5бестей! [ауег; 4 - пеа-т;зшайпе [ауег, 
5 - Пей! сагтег осай2аноп агеа 
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Рис. 2. Область расчета цилиндрического сегмента коллектора, 
в случае, если источник света расположен под прямым 
углом к диаметру сегмента коллектора 


Ето. 2. Сасшайоп аотат офсуйпатса[ соПесюг 5естеп! т са5е 
уйеп Пой! зоигсе 15 потта[ ю Фатеег оф соПесюг бестеп! 


Математическая модель задачи имеет вид: 


область О! а (2 этааТ °) 0 (1) 
область 2› ав (^, этааГ °) =0. (2) 
Примем следующие граничные условия 
р К, › 
эт?) 
2 — о |7 (&,ф)-Т, |; 
р=^,, 
эт эт) 
7 (®,,ф,!) — 72) (®,,ф,!) 5 2 бр. — р бр. 5 


р = ^.. 
Для расчета температурного поля будем использовать метод конечных элементов (МКЭ). На рис. 3 
представлена конечно-элементная сетка, покрывающая область расчета. 





Рис. 3. Конечно-элементная сетка области расчета 


Ето. 3. Етйе-@етет! ета ор сасшапоп аотат 
Решение данной задачи сводится к итерационному численному решению прямых задач и минимизации 
функционала 
м 
=> (т,-Т“"), (3) 


т Л у] 


где М — число точек измерения температуры Т;”, Т; — значения температуры в тех же точках, полученных 
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путем решения прямых задач методом конечных элементов. 

МКЭ в настоящее время стал одним из самых распространенных и эффективных методов решения 
задач математической физики и техники, сформулированных в виде краевых задач для дифференциальных 
уравнений с частными производными. 

Пусть требуется найти распределение температуры Т(№) в некоторой области Г, удовлетворяющей 
уравнению стационарной теплопроводности 


а» зтаат(\)} + ГСМ) =0 (4) 

и условию на поверхности 5, ограничивающей область Г, 
»5-(м) =а(м)т(м). (5) 
где (№ ) — коэффициент теплопроводности, в общем случае функция точки МЕЙ; о (М ) — коэффициент 


теплоотдачи, Мех; Г (№ ) — объемная плотность источников тепла. 


Первый этап МКЭ — построение неортогональной сетки, покрывающей область расчета Г. Сетка 
образует совокупность конечных элементов, например, тетраэдров, в плоскопараллельном поле треугольников. 
Далее следует нумерация элементов и узлов, составление массивов координат узлов. 

На втором этапе выполняется аппроксимация искомых и известных функций (в нашем случае Т(М) и 


Уд (№ ) ). На третьем этапе будет показано, что математическую модель исходной задачи можно преобразовать 


так, что аппроксимацию функций Т(№и / (№ ) допустимо искать в пространстве дифференцируемых функций 


п 
С. Построим в этом пространстве подпространство с конечным базисом К Базисные функции \у, (№ ) 
выберем финитными (с конечным носителем) и кусочно-линейными. 
В/-м узле сетки положим \у, =1, в остальных узлах \у, =0. 


Тогда аппроксимирующие функции на втором этапе МКЭ будут иметь вид 
Г(жу,2)= ХТ, (х,у,2); У (х,у,2)= > Дм, (х,у,2), 


где и — количество узлов сетки; 7’, = Г 2. ‚5 Л, = а й 


В пределах А-го конечного элемента искомая функция представляется в виде 
> т-4 В 
Г” (х,у,2)= > Тру, (х,,2), (6) 
=т-1 


К р р р 
где \у, =а;х+Б, ус, 2+а,. 
В случае плоскопараллельного поля имеем 
т — т-3 | 
Г" (х,у)= > Ту, (х,У), 
1=т-1 
К К 
где лу, = агх+Ь:у+с, 
[ч 
Постоянные коэффициенты \у, определяются через координаты узлов элемента. 


На третьем этапе интегралы по области 7 заменяются суммой интегралов по конечным элементам И”, 
интегралы по 5 — суммой интегралов по граням элементов, аппроксимирующих ъ. 

Строится система уравнений МКЭ и учитываются граничные условия. Эту систему называют 
дискретной математической моделью задачи (ДММ). 

Рассмотрим построение ДММ для проекционного варианта МКЭ, суть которого заключается в 
минимизации нормы невязки 


== Ал\етааТ (№) + Л 
И 
путем ортогонального проектирования ее на конечномерное подпространство с базисом | ‚. Минимум 
= 


невязки обеспечен в том случае, если все скалярные произведения невязки # и \у, будут равны нулю. 


В МКЭ принято использовать скалярные произведения более широкого класса функций ([>). 
В результате получим систему МКЭ в виде 
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у |4’ (^етааТ)+ У ай =0 ‚1=1,2....П. 


Используя известные тождества 


(фа) = фАгла + аетааф ; 
Паг»(фа)аТ = [фа„ 45, 


где ф — скалярная функция; а — векторная функция; а, — проекция а на внешнюю по отношению к Г 


нормаль, получим систему уравнений вида 
- 7 
Гл ета4Тэгаауу АУ = [[лу, тазу, У Е не (7) 
| 5 И | 


где и — число узлов сетки, в которых неизвестны Г; 


Дальнейшие преобразования (7) приводят к системе уравнений вида 
2% Вт, = 2.4 [Гугаб+ хх чичаТ , ый (8) 
г. ты Я |7 


где Ч = р 
А 


— нормальная производная Г (М ) на грани 5 Еб 


> К К К 
А-го конечного элемента (принимается постоянной); И” — область А-го элемента; В, = | эгай\у; зта4у ‚АТ. 
у* 


Внешнее суммирование в (8) (по А) ведется по всем элементам, имеющим общий узел 1, внутреннее (по 
7) — по всем узлам А-го элемента. Для всех внутренних узлов первый интеграл в правой части (8) равен нулю. 

Систему (8) называют дискретной математической моделью задачи, если в нее подставить граничные 
условия. Решение системы (8) и использование (6) позволяют определить все интересующие величины 
(градиенты, потоки и т.п.). 


Исходные данные: температура окружающего воздуха Т,,=23 °С; тепловой поток 4=4, $ШФ, 
4, =319 Вт/м” (значение теплового потока от источника излучения). Теплопроводность: светопропускающего 
слоя Л =0,03 Вт/мК;  теплопоглощающего слоя Л, =0,34 Вт/м-К; коэффициент теплоотдачи 
о... =9,5 Вт/м”.К; К =0,05 м; К, =0,0492 м; К, =0,0442 м; ф =30`; ф, =60`; ф. =90'. 

Необходимо определить распределение температур в установившемся режиме в точках: Т'(К,,ф,); 
Т,(®,ф.); 7, (К.Ф); Т,(В,Ф,); 1, (К.Ф, ); Ту (В,Ф,). 


Результаты вычислительного эксперимента приведены в таблице] : 


Таблица 1 
Табе 1 


2%) | (8%) | 58%) | тв) | т) | 18%) 


Вычисленные 
48,7 41,3 56,7 46,8 60,1 49,2 





значения, °С 


Далее приведены результаты математического моделирования теплообмена в одном сегменте 
коллектора для случая, когда источник света расположен под углом 45° к диаметру образца (рис. 4). 
Исходные данные соответствуют тем, которые использовались при моделировании первого случая. 


Найдем распределение температур в установившемся режиме в точках: Т(К,,ф,); Т,(К,ф,); Т, (В,,ф.); 


Т, (А,,$,); Т. (В, $); Т, (В, ); Т, (В,,ф. ); т, (®,ф.); Т,(В,ф5); По выде 
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Рис. 4. Область расчета для случая, когда поток солнечного излучения расположен под углом 
к диаметру сегмента коллектора 
Ето. 4. Сасшаноп аотат юг сабе ийеп 5о[аг гайапоп Пих 15 5её а! 45° апёе 
0 соПесюг 5есотет! йатеег 
Значения плотности теплового потока, приведенные в таблице 2, найдены на основе решения обратной 
задачи теплообмена. 

Таблица 2 

Табе 2 


1020 20-30 | 3060 | 60—70 70-80 [80-90 | 90-100 |100—110|110-120|120-130| 130—180 
|182 | 194 | 202 | 210 | 202 | 194 | 182 | 126 | 12 | 17 | Ш! | 105 
Вт/м 


Адекватность математической модели была проверена экспериментально. С этой целью был 





изготовлен замкнутый цилиндрический сегмент коллектора. В качестве абсорбера использовалось шерстяное 
полотно. В ходе эксперимента образец воздушного коллектора был установлен вертикально. Напротив образца 
на расстоянии 0,8 м располагался источник света — прожектор с галогеновой лампой мощностью 1кВТт. 
Излучение на поверхность теплопоглощающего слоя было направлено под углом 45°. С помощью датчиков 
температуры ИВТ (прибора для измерения плотности теплового потока) были зафиксированы значения 
температуры наружной поверхности светопропускающего и теплопоглощающего слоев в период 
термостабилизации, т.е. по истечении 20 мин. воздействия на коллектор источника света. 
Результаты измерительного и вычислительного экспериментов приведены в таблице 3. 
Таблица 3 
ТаЫе 3 


р 5) | 55) [ а) | тв) 
Измеренные 477,0 39,0 48,0 36,0 35,6 
значения, °С 


Вычисленные 46,1 38,2 44,0 36,7 38,8 


значения, °С 
р | К | пех | 1 | 5 | 


Измеренные 30,6 33,0 30,0 33,4 ьь 
значения, °С 


Вычисленные 33,1 30,0 


значения, °С 


-82 т 
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Следовательно, результаты моделирования процессов теплообмена в сегменте коллектора для разных 
направлений воздействия потока солнечного излучения согласуются с результатами измерений с достаточной 
для инженерных задач точностью. 
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Выводы. Прирост температуры теплопоглощающего слоя по отношению к температуре окружающей 
среды составил от 7° до 26,2° в зависимости от рассматриваемой части сегмента коллектора и направления 
воздействия потока солнечной радиации. Таким образом, локализуясь между теплопоглощающим слоем и 
теплоизоляционной подложкой, конструкция коллектора обеспечит нагрев теплоносителя. Максимальные 
значения температурного поля поверхности теплопоглощающего слоя (60°) коллектора отмечаются при прямом 
воздействии солнечного излучения. Изменение угла воздействия потока излучения до 45° приводит к 
снижению температуры. В зависимости от сектора воздействия, температура поверхности теплопоглощающего 
слоя может снижаться на 3—14° по отношению к замерам при прямом воздействии солнечных лучей. 

Адекватность разработанной математической модели позволит разработать программный продукт, 
который автоматизирует процесс расчета температурного поля коллектора. 
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